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The discovery of many novel realizations of the inflationary universe paradigm has led to a de- 
generacy problem: many different inflationary Lagrangians generate the same perturbation spectra. 
Resolving this problem requires the future discovery of additional observables, beyond the scalar 
adiabatic and tensor two-point functions on CMB scales. One important source of degeneracy arises 
in models where the density perturbation is generated by a non-inflationary degree of freedom, for 
example, through curvatons or modulated reheating. We consider the curvaton scenario as repre- 
sentative of this class, and analyze the degeneracy with single field, canonical inflation that results if 
the curvaton goes undetected by future observations. We perform Monte Carlo potential reconstruc- 
tions in the absence of distinguishing observables, such as non-Gaussiantities or isocurvature modes. 
The resulting degeneracy is considerable and the improved measurements of spectral parameters 
from future probes like CMBPol, offer little to better the situation. Given a degeneracy-breaking 
observation, the observables must still be inverted to obtain the inflationary potential, with differ- 
ent observations resulting in reconstructions of varying quality. We find that a future detection of 
isocurvature modes or a precision measurement of the tensor spectral index will enable the most 
successful reconstructions in the presence of curvatons. 



I. INTRODUCTION 

Observational cosmology has determined that the uni- 
verse is homogenous and isotropic on scales larger than 
100 Mpc and that the universe is expanding according 
to Hubble's law. Measurements of the cosmic microwave 
background radiation (CMB) reveal that this homogene- 
ity and isotropy existed even at the time of recombina- 
tion to an extreme accuracy at the order of 10 -5 on all 
scales up to the present horizon. The modern standard 
model of cosmology attributes these observed properties 
of the universe to an early accelerating period, known as 
inflation [1-3]. 

The inflaton, the field that dominates the energy den- 
sity of the universe during inflation, traditionally as- 
sumes a dual role: it drives the accelerated expansion 
and generates the primordial fluctuations that seed large 
scale structure. Precision measurements of the CMB 
and large scale structure surveys (LSS) have been in- 
strumental in constraining the form of the inflationary 
Lagrangian. This process of utilizing cosmological obser- 
vations to limit the form of the inflationary Lagrangian is 
called reconstruction, and relies on the assumption that 
there exists a unique mapping between the set of observ- 
ables and the free parameters of the Lagrangian. Such a 
mapping exists in the framework of single field, canonical 
inflation under the slow roll approximation, and the pro- 



gram of reconstruction has been extensively developed 
[4-7]. Of course, this process assumes that the inflaton 
is responsible for generating the perturbations; if this 
assumption is invalid then any process that seeks to 'in- 
vert' a subset of observables to obtain the underlying 
free parameters of the Lagrangian will reveal a space of 
Lagrangians that is observationally degenerate 1 : many 
distinct theories will provide the same subset of observ- 
ables [9]. 

It is possible to liberate the inflaton from its role in gen- 
erating the primordial perturbations. In curvaton scenar- 
ios [10-12], this role is assumed by an additional light de- 
gree of freedom decoupled from the inflationary dynam- 
ics; in modulated reheating, perturbations are generated 
by fluctuations in the decay rate of the inflaton [13, 14]. 
In these theories, the central assumption of traditional 
reconstruction-that the inflaton generates the primordial 
spectra-is violated. Without knowledge of how the spec- 
tra were generated, whether by the inflaton or by some 
other means, a unique inversion of observables is clearly 
impossible. The key to resolving this impasse is to recog- 
nize that while degeneracies might exist within a subset 
of observables, it might be possible to utilize additional 
observables as a means of distinguishing between differ- 
ent models. While this is the obvious solution, it is not 
clear which observables are relevant to which scenarios, 
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We refer here to the degeneracies in the locally reconstructed 
inflaton potential; degeneracies of a different type arise globally 
as discussed in [8]. 
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or how effectively such observations enable reconstruc- 
tion. It is now possible to state the degeneracy problem: 
what observations must be made to distinguish between 
different theories, and if these observations are not made 
how large is the resulting degeneracy? The subsequent is- 
sue, one that we spend considerable time investigating in 
this work, is the inversion problem: given a distinguish- 
ing observation, how well can the set of observables be 
inverted to obtain the underlying inflationary potential, 



II. POTENTIAL RECONSTRUCTION IN 
CANONICAL INFLATION 

The program of potential reconstruction was initiated 
in 1990 [4], and subsequently expanded and refined [5-7]. 
The reconstruction program initially targeted inflation 
driven by a single, canonically normalized scalar field, 0, 
minimally coupled to gravity, 



(1) 



In this paper, we examine these questions within the 
context of the curvaton scenario as a prototype of the 
case in which the perturbation spectra are not generated 
by the inflaton. We do not consider modulated reheating, 
as it is phenomenologically similar to the curvaton. Our 
analysis employs Monte Carlo reconstruction, in which 
the inflationary model space is stochastically sampled 
and models satisfying specific observational constraints 
can be identified and extracted for further analysis. We 
perform reconstructions for single field canonical inflation 
(serving as a reference model) and the curvaton scenario. 
We first perform the reconstruction in the absence of any 
distinguishing observables to determine the size of the de- 
generacy, and then systematically investigate the effects 
that an observation of isocurvature perturbations, non- 
Gaussianities, and a precision measurement of the tensor 
spectral index each has on reducing the degeneracy and 
enabling reconstruction. 



This analysis is far from exhaustive; there are scores 
of inflationary scenarios beyond single field inflation and 
curvatons that could conceivably contribute to the de- 
generacy problem. For example, the altered dynamical 
degrees of freedom of non-canonical inflation [15] also 
contribute to the degeneracy problem [16], and we inves- 
tigate Dirac-Born-Infeld (DBI) inflation [17, 18] in detail 
as an example in [19]. While other sources of degener- 
acy exist, we consider the curvaton and DBI inflation to 
be two well-motivated, well-studied and phenomenologi- 
cally rich scenarios that are representative of two strong 
sources of degeneracy: non-inflaton produced primordial 
spectra and non-canonical inflationary dynamics, respec- 
tively. In the absence of a fully comprehensive analysis, 
the results presented here combined with [19] represent 
a 'best-case' scenario for reconstruction, and should be 
viewed as a first step towards understanding how large 
the degeneracy might be, and which cosmological observ- 
ables are needed to mitigate the problem. Despite the 
significant challenges we identify in this study, it is our 
goal to determine which feasible cosmological observa- 
tions will help us to best understand the fundamental 
physics behind the inflationary universe paradigm. The 
information gathered in this analysis is crucial to the 
planning and optimization of future cosmological exper- 
iments. 



The analysis provides a framework for mapping power 
spectrum observables to the coefficients of a Taylor ex- 
panded inflaton potential. The method begins with the 
Hamilton- Jacobi equation, 
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where H(cj)) = a /a is the Hubble parameter, a(t) is the 
scale factor and / = d/d(j). Next, the Hubble parameter 
is Taylor expanded about some field value, (/>o, 
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At zeroth-order, H((j)) is a constant, corresponding to the 
de Sitter solution. The first, second, and higher order 
terms are [22] 
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where the functions e = 2Mp l (H' / H) 2 , rj = Ai, and A n 
(n > 2) parameterize deviations from pure de Sitter ex- 
pansion. These are commonly referred to as slow roll 
parameters in the literature; however, they are defined 
here without any assumption of slow roll and will hence- 
forth be referred to as flow parameters. Combining Eqs. 
(4), (3), and (2) and collecting like-powers of (0 — 0o) 
leads to the Taylor coefficients of the potential, 
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We note that these expressions are exact to all orders in 
slow roll. 

Potential reconstruction reveals the fact that devia- 
tions from pure de Sitter are manifested in the power 
spectrum of density fluctuations as deviations from scale 
invariance. The power spectrum of curvature perturba- 
tions can be parameterized as a power law, 
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P*(k) = P*(ko) [ — 



(9) 



where k is the comoving wavenumber of the fluctuation 
(with fco defined as the scale that crosses the horizon 
when <j) — 0o ) an d n s is the spectral index. Deviations 
from pure de Sitter also result in the generation of a 
large-scale gravitational wave, or tensor, spectrum: 



P h (k)=P h (k ) 



k 



(10) 



To first-order in slow roll, the amplitude of the power 
spectra at horizon crossing (k = aH) are given by the 
well-known expressions, 
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The scale dependence of the spectra are then determined 
by the time variation of H((j)) and e(0): 

_ dlnP^(k) 



n s (k ) 



n T (k ) 



dink 



d\nP h (k) 
dink 



2r](M -4e(0o), 



-2e(0 o ). 



(13) 
(14) 



At first order in slow roll, the flow parameters e and rj 
enable one to map the spectrum observables, Eqs. (11- 
14), to the Taylor coefficients of the inflaton potential, 
Eqs. (6-8). For example, at lowest order, we have 
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(15) 
(16) 

(17) 



where r = Ph(ko)/P$(ko) is the tensor /scalar ratio, and 
we have neglected all higher-order terms in Eqs. (6 - 
8). While only first-order approximations, Eqs. (15 -17) 
provide a wealth of information. For example, given the 
current measurement of the amplitude of scalar pertur- 
bations, P$(/uq), a detection of tensors by ESA's Planck 
Surveyor [23] will determine V(4 > o) an d V'^o) to within 
a percentage error of roughly Ar/r. The prospect of re- 
constructing V((j)o) is especially exciting because it cor- 
responds to the energy scale of inflation - an important 
clue to uncovering the identity of the inflaton. 
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FIG. 1. The 'zoology' of inflation models. 



The reconstruction equations (15 -17) reveal that in- 
flationary potentials can be grouped into three distinct 
classes based on their observable predictions (c.f. Fig- 
ure 1). The economy afforded by this classification - 
that vast numbers of different inflationary potentials or- 
ganize themselves into a few observational families - 
came to be known as the 'zoology' of inflationary mod- 
els [24, 25]. Models labeled 'hybrid' include potentials 
that evolve asymptotically to their minima, requiring an 
auxiliary field to end inflation. However, they are ef- 
fectively single field models with non-vanishing energy 
density at the minimum, and have the general form 
V((j)) oc 1 + ((/>//i) p , where \i is an energy scale and p 
a positive integer. Hybrid models are characterized by 
the conditions V"(<j>) > and (logV(0)) // > 0; the sim- 
plest models of tree- level hybrid inflation [26, 27] belong 
to this class. 'Small field' and 'large field' models are dif- 
ferentiated by their initial field values. Large field models 
are characterized by a field initially displaced far from its 
minimum, with the general form V(<j)) oc (^)//i) p , satisfy- 
ing V"(<j>) > and (logV(0)) // < 0; m 2 (j) 2 chaotic infla- 
tion [28] is one example. Conversely, small field models 
are characterized by a field initially close to the origin, 
with general form V{4>) oc 1 — (0//i) p (near the maxi- 
mum), satisfying V"{(j)) < and (logV(0)) // < 0; 'new' 
inflation and other models based on spontaneous sym- 
metry breaking belong to this class. While phenomeno- 
logically distinct, these classifications are also useful for 
discerning more fundamental aspects of inflation. For 
example, large field models support chaotic initial condi- 
tions and eternal inflation, but are difficult to embed in 
a string theoretic framework; 2 small field models with 
A(j) < Mpi are well-behaved effective potentials, but must 
have fine-tuned initial conditions; hybrid models are mul- 



See, however, [29, 30] 
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tifield scenarios. 

While future observations promise to reduce the size 
of errors in the n s -r plane, it will still be possible to 
generate ensembles of potentials that provide good fits 
to the data. 3 However, with any luck, we will be able 
to determine which class of potential, of the types dis- 
cussed above, is ultimately responsible for driving infla- 
tion. This of course all assumes that the simplest imple- 
mentation of single field, canonical, minimally coupled 
inflation is true. To what degree is our ability to recon- 
struct the physics of inflation threatened by relaxing one 
or more of these assumptions? In the next section, we 
analyze the effects of a curvaton on the reconstruction 
program. 



III. CASE STUDY: RECONSTRUCTION IN 
THE PRESENCE OF CURVATONS 

A. The Scenario and the Power Spectrum 

The curvaton mechanism [10-12] relaxes the assump- 
tion that the initial curvature perturbation originated 
during inflation, and instead attributes it to the fluctu- 
ations of a second field, a, known as the curvaton. The 
curvaton is a weakly coupled scalar field that is relatively 
light during inflation, m 2 <C H 2 . After inflation ends, the 
initially displaced curvaton rolls to its minimum where 
it begins to oscillate during the post-inflationary radia- 
tion dominated phase. These oscillations set up a small 
isocurvature perturbation that grows as the curvaton 
continues to oscillate. After the curvaton decays, the per- 
turbation is converted to an adiabatic mode and struc- 
ture begins to evolve according to the standard model. 
The curvaton model is phenomenologically rich; depend- 
ing on the energy density of the curvaton at the time of 
decay, there may be residual isocurvature modes or pri- 
mordial non-Gaussianity large enough to be detected by 
future experiments. 

We first discuss the effect of the curvaton on the pri- 
mordial power spectrum. In general, both the curvaton 
and the inflaton will contribute to the final curvature 
perturbation, 
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where $ is the Bardeen potential (curvature perturba- 
tion in longitudinal gauge), Scj) and 5 a are the inflaton 
and curvaton vacuum fluctuations, and denotes the 
value at the end of inflation. The function f(a) parame- 
terizes the amount that the curvaton contributes to the 
curvature perturbation. 

In what follows, we will drop the subscripts and 
with the understanding that all spectral observables 



are given during the radiation dominated phase follow- 
ing curvaton decay. The curvaton is effectively massless 
during inflation so that, 5 a = H/2ir, and the curvature 
perturbation Eq. (18) becomes 
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Since the inflaton and curvaton perturbations are uncor- 
rected, one can write the power spectrum [35-37], 



1 + 
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where / = Sf/y/2. The spectral index follows, 

1 <flnP*(fc) 2r?-2e 

n s — 1 = =— ; — = — 2e H = . 
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(20) 



(21) 



Because the energy density of the curvaton is strongly 
subdominant during inflation it does not generate grav- 
itational waves, and the tensor spectrum is the same as 
that generated in single field inflation. However, the ten- 
sor/scalar ratio is modified as a result of the new scalar 
amplitude Eq.(20), 
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Note that by taking / — >> we recover the usual inflaton- 
generated spectrum Eqs. (11) and (13), and so by tuning 
/ we can control the contribution of the curvaton to the 
final perturbation. The 'pure curvaton' limit is achieved 
when f 2 (cr)e — >> oo. 

The addition of the curvaton alters the relationship 
between the observables r and n s , allowing for a richer 
phenomenology. By relegating responsibility for generat- 
ing a nearly scale invariant power spectrum to the curva- 
ton, constraints on the form of the inflaton potential are 
relaxed [37-39]. However, it is also apparent that the in- 
clusion of the curvaton threatens a unique reconstruction 
of the inflaton potential from observations. The observ- 
ables r and n s pick up a dependence on /(cr), and pass 
it along to the potential coefficients. Solving for e and rj 
in Eqs. (21)-(22) and using these relations in Eqs. (6-8), 
we obtain new expressions for the potential coefficients, 
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See, for example, [31-33]. 



Comparing Eqs. (15-17), we see that V((/)q) is unchanged 
by the inclusion of the curvaton because its energy den- 
sity is subdominant to that of the inflaton. However, 
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there is no longer a unique inversion from the observables 
P$(fco), r i an d n s to the coefficients V(<po), V'((/)q), and 
V"((j)o), and the degree of resulting uncertainty depends 
on the possible values of f(cr). In the limit f 2 (cr)e ^> 1, 
the curvaton contribution dominates the overall curva- 
ture perturbation Eq. (19), with the result that 

r^J^, (26) 

n s - 1 ->> -2e. 
In this limit, the potential coefficients become 

7T 2 



V(^o) = yM;Pi(l„)r(i ), 

3 V(0o) f dn s /d\nk(ko) 



(27) 

(28) 
(29) 



t 7-// / / \ 3V(^) fdn s /dlnk(k ) \ 

where dn s /d\nk is the spectral index running. A com- 
parison between these expressions and Eqs. (15-17) re- 
veals that by tuning f(cr) from 'no curvaton contribution' 
(single field) to 'maximal contribution' (pure curvaton), 
we obtain very different relations between the inflaton 
potential and the observable parameters, with the possi- 
bility of the same potential giving rise to a wide range 
of observables, depending on the value of f(cr). This is 
the degeneracy problem that threatens reconstruction in 
the presence of curvatons. Equations (23-25) are derived 
from the lowest order expressions for the observables r 
and n s ; we neglected terms of order e 2 , eq, etc. in ob- 
taining Eqs. (21) and (22). These expressions thus serve 
merely to illustrate the emergence of possible degenera- 
cies between single field and curvaton models; we conduct 
a robust, higher-order reconstruction in § HIE. 

In this analysis we allow a substantial range of val- 
ues for /(a), including the limit in which the curvaton 
contribution is large. This limit is attained if the cur- 
vaton begins to oscillate early relative to the time that 
it decays; this way, the initially small curvaton fluctua- 
tions have time to grow. Oscillations must then begin 
when the curvaton energy density is strongly subdomi- 
nant: H 2 ~ m 2 ^> 4j2-<J 2 , and therefore the initial field 
value must satisfy a <C Mp\. In this case, under the 
sudden decay approximation [35], 

K°) = rA — , (<x«M P1 ), (31) 
y cr 



where 
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and ^cr,dec is the density parameter of the curvaton, 
Q a = at the time of decay 4 . The contribution 



4 A more realistic gradual decay of the curvaton leads to similar 
results [40, 41]. 



to the overall curvature perturbation is maximized when 
the curvaton comes to dominate the universe before de- 
cay, in which case r a = Q a ,dec = 1- 

Note, that for a ~ Mpi, the curvaton oscillates late; 
the contributions of the curvaton and inflaton to the fi- 
nal curvature perturbation become comparable, and f(a) 
must be obtained numerically [36, 42]. However, even for 
a <C Mpi, the curvaton can be subdominant if r a <C 1. It 
is also possible that a ^> Mpi, leading to a second bout of 
curvaton-driven inflation after the inflaton energy density 
drops sufficiently. This case is essentially double inflation 
[43, 44], and we will not discuss it further. 

The curvaton scenario is phenomenologically rich, with 
the ability to generate large non-Gaussianity and a sig- 
nificant isocurvature perturbation. A detection of either 
of these would serve to break the degeneracy with single 
field inflation. The problem then becomes one of inver- 
sion: if we are unable to constrain the function /(cr), then 
potential reconstruction will be degraded. 



B. Isocurvature Modes 

As a multifield model, the curvaton scenario is capa- 
ble of producing an isocurvature perturbation in addition 
to the adiabatic perturbations discussed in the previous 
section. A generic isocurvature perturbation between two 
separate fluid components with energy densities pi and 
Pj can be written, 



(33) 



where we have introduced the curvature perturbation, £, 
which is given during matter domination by ( = — 5<E>/3. 
The curvaton can generate isocurvature perturbations in 
the decoupled fluids after inflation: cold dark matter 
(CDM), baryons, and neutrinos [35], 
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where 7 denotes radiation. We define the power spectra 

(27T) 3 S(k 1 +k 2 )P C (k 1 ), 

(2n) 3 5(k 1 +k 2 )P s (k 1 ), 
(27r) 3 5(k 1 +k 2 )P c<s (fc 1 ). 
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Most cosmological observations are sensitive to the rel- 
ative amplitudes of the isocurvature and adiabatic spec- 
tra, |<S/C| 2 , as well as the degree of correlation between 
the two modes. The former quantity is conventionally 
measured in terms of the parameter, 
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A standard CDM iso curvature mode can be generated 
in either of two ways, depending on whether the CDM 
was created before curvaton decay, or whether it was cre- 
ated as a direct product of this decay. In the former 
case, the CDM perturbation is generated by the inflaton, 
Ccdm = Cinf, and we obtain 



£cdm 
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where A = / 2 (<j)e, we have taken £ 7 = £, and Eq. (20) 
has been used to relate the inflaton fluctuation to the 
overall curvature perturbation. In the pure curvaton 
limit (f 2 (a)e — >> oo and thus A ^> 1), the isocurvature 
contribution is too large and has long since been ruled 
out [35]. However, a sufficiently large inflaton contribu- 
tion can reduce the amplitude of <Scdm an d bring it into 
agreement with current bounds; this type of isocurvature 
perturbation is uncorrelated with the adiabatic compo- 
nent ((3 = 0), with the constraint a < 0.077 (95% CL) 
[46]. If future CMB missions observe a CDM isocurva- 
ture perturbation then single field inflation will be ruled 
out and the degeneracy broken, and then a detection of 
tensors will suffice to enable a successful inversion 6 - we 
will have four functions (V(0o), V'((/)q), V"((j)o), f(cr)) in 
terms of four observables (P$(fco), r, n s , and ao). The 
Planck surveyor is expected to obtain a percentage error 
of around 4% for CDM perturbations at 68% CL [47], and 
we will analyze this case in detail in the next section. 
However, if the CDM is not generated before curvaton 
decay, then no such isocurvature mode will be produced. 

The other possibility is that the CDM is created as a 
direct product of curvaton decay. In this case [48] 
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(43) 
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where r a is as defined in Eq. (32). In the pure curvaton 
limit (A -1 — » 0), the isocurvature perturbation is anti- 
correlated (/? = —1), and its amplitude has been con- 
strained in several analyses [49-54], most recently using 



WMAP7+BAO+SN [46] giving a_i < 0.0037 (95% CL). 
From Eq. (22), it is possible to relate the tensor/scalar 
ratio to the correlation angle, 
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and so anti-correlated CDM isocurvature modes corre- 
spond to a lack of tensors, r = 0. The loss of r as an 
observable precludes a determination of the energy scale 
of inflation, V(<j>o)', however, the terms V /V and V" /V 
are completely determined by n s and dn s /dlnk (c.f. Eqs. 
(28-30)). We will investigate this case in the next section. 

While the un- and anti-correlated isocurvature pertur- 
bations discussed above are given the most attention in 
the literature, it is evident from Eq. (43) that if the cur- 
vaton and inflaton contributions to the curvature pertur- 
bation are both relevant, than it is possible to obtain a 
CDM isocurvature mode with an arbitrary correlation 7 . 
In this event, given measurements of the observables: r, 
n s , and P$(fco), a constraint on j3 will determine f(a) 
and enable a successful reconstruction, with Planck ex- 
pected to achieve A/3 ~ 0.04 (65% CL) [47]. We will 
investigate this possibility in the next section. 

Before concluding this subsection, we mention that 
baryon isocurvature fluctuations appear identical to 
CDM modes in the CMB, apart from a factor of 
Qb/QcDM- Therefore, the analyses that we carry out 
in §111 E for the above cases involving CDM isocurvature 
perturbations apply to baryonic modes as well. 



C. Non-Gaussianity 

In addition to isocurvature modes, the multifield set- 
ting of the curvaton scenario has the potential to gener- 
ate large non-Gaussian temperature fluctuations in the 
CMB. The non-Gaussianity of the fluctuations give rise 
to a nonzero bispectrum of the curvature perturbation, 



(CCO^(2^(ki + k 2 



k 3 )^ 0t (fcl,fc2,fc3), 



(46) 

where C tot denotes the total curvature perturbation in the 
matter dominated era, and comprises both adiabatic and 
isocurvature components. The bispectrum thus contains 
the four terms: (CkiCk 2 Ck 3 )> (CkiCk 2 5 k3 >, (Cki<S k2 S k3 ), 
and (<Ski<$k 2 <Sk 3 ) [48, 56, 57]. In order to make contact 
with the parameter f^L that is constrained in contem- 
porary analyses, we define 



Our sign convention is such that the temperature fluctuation 
AT/T = C/5 - 2<S/5 (take C = K -K = K to compare 
with WMAP [45]). Physically, a correlated perturbation ((3 = 1) 
corresponds to a suppression of the temperature power spectrum 
at the low multipoles. 

This assumes that the curvaton can be distinguished from other 
mechanisms that might generate a CDM isocurvature mode, such 
as the axion. This might require the simultaneous detection of 
non-Gaussianities consistent with the thermal scenario consid- 
ered. 



k 2 , k 3 ) = -/jvl [Pc{ki)P^k 2 ) + 2 perms] , 



(47) 



7 See [55] for the most recent constraints on arbitrarily correlated 
isocurvature perturbations. Their sign convention on (3 is oppo- 
site ours. 
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with y NL = l(f£l + f$p « If^r = 6 CCC + |6 CC5 + 

ij^CSS + ^y^sss where fffi L is related to the three-point 
function of the fields at horizon crossing, and /jyff 1 
arises from the nonlinear relation between the curva- 
ture and field perturbations; as we have indicated, the 
'local' form dominates the bispectrum in the curvaton 
scenario. The nonlinearity parameter, /nl, character- 
izes the deviation of the fluctuation from pure Gaussian, 
C = Cg- \Inl {Cg ~ ((g)) j where ( g is the Gaussian cur- 
vature perturbation. Since it is local in real space, the 
Fourier modes satisfy ki « k2 ^> IC3. Using second-order 
perturbation theory under the sudden decay approxima- 
tion, the non-Gaussianity arising from the adiabatic com- 
ponent is [58-60], 



^(adi) 5 



5 1 (| - 2r a 



(48) 



The non-Gaussianity arising from the isocurvature com- 
ponent, fj^L^ — \yj^sss^ and that of the cross- 
correlated components oc b^s and b^ss, depends on the 
thermal history of the universe: for the different CDM 
isocurvature modes considered in the previous subsec- 
tion we have bff L s = -3b% s £ = %% c £ = -276$£ for 

~„°NL 



the uncorrelated mode and bffjf 
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r^^NL ^ or ^ ne anti-correlated mode. Evidently, a suffi- 
ciently low curvaton energy density at the time of de- 
cay (r a <C 1) will result in large and potentially ob- 
servable non-Gaussianities, and it is therefore possible 
for the non-Gaussianity produced by the isocurvature 
mode to dominate the signal. However, in the pure cur- 
vaton limit it is in general difficult to generate large 
non-Gaussianities without simultaneously violating the 
bounds on the isocurvature amplitudes. For example, in 
the pure curvaton limit, the bound on the anti-correlated 
mode (a_i < 0.0037) corresponds to r a > 0.98 and a 
bound G [-1.25,0.43]. This lies below the expected 

limits of future missions: the Planck surveyor is expected 
to reach a 68% limit of A/]^£ al = 5 [61] while large scale 
structure surveys should yield A/jy£ al ~ 1 — 5 [62, 63] 
at 68% CL. Simultaneously observable non-Gaussianities 
and isocurvature contributions are, however, possible 
for certain neutrino [64] and baryonic modes in the 
case of mixed curvaton and inflaton perturbations [65]. 
However, as discussed in the last subsection, a detec- 
tion of isocurvature modes independently imposes a con- 
straint on the function /(a), and measurements of non- 
Gaussianities are not needed for reconstruction. A detec- 
tion of /nl is, however, especially advantageous in the 
event that no isocurvature perturbations are produced in 
the scenario because it will break degeneracies and might 
help with inversion. 

In the absence of isocurvature production, the only 
nonzero contribution to the non-Gaussianity is from 
/jvx^, which is a function of both r G and f(cr). In or- 
der to constrain /(a), an additional, complementary ob- 
servation is needed. In addition to the bispectrum, the 



non- Gaussian fluctuations generate a trispectrum statis- 
tic [66-68], 

(C kl Ck 2 Ck 3 Ck 4 ) = (2k) 3 5 Td^iMMM), 

(49) 

where, in contrast to Eq. (46), we include only the adi- 
abatic perturbations since we are considering the case in 
which isocurvature modes are absent. 8 The trispectrum 
can be written in terms of the estimators tnl and g^L 
[70, 71], 

T c (fci, fc 2 , &3, h) = TNLiP^P^k^P^h) + 11 perms) 
54 

7gN L (Pdk 2 )P c (k 3 )P c (k 4 ) + 3 perms). 

(50) 



25 - 



The function g^L contains higher-derivative terms and 
is subdominant relative to tnl [41], and we will not con- 
sider it further. The estimator tnl can be related to the 
physical parameters of the model using the SN formalism 
[68, 72, 73], 



36 
25 J 



TNL = ^f NL (l + X- 1 ) 2 . 



(51) 



Happily, the trispectrum does not introduce additional 
free parameters, and a measurement of f 1 ^ 1 and r 1 ^ 1 
will constrain the functions f(cr) and r a . The 95% limit 
on the bispectrum from WMAP7 is -10 < < 74 

[46], while the trispectrum is much less well-constrained: 
\tnl\ < 10 8 from COBE [74]. Future full-sky galaxy 
surveys might achieve an accuracy of A/]^£ al ~ 1 — 5 
(68% CL) [62, 63] and Planck is expected to achieve 
At n c l 1 ~ 560 ( 95 % CL ) [71]. We will study the effect of 
such detections on reconstruction in §111 E. 

In general, it is possible for the amplitude of non- 
Gaussianity, /jvL a1 ' t° vary with scale, /nl ~ k nNG , 
where the spectral index is defined 



n N G 



NL 



dink 



(52) 



Scale dependent local non-Gaussianities [75] can arise in 
the curvaton scenario if the curvaton potential has inter- 
action terms [76] or if the final density perturbation is 
a comparable mixture of inflaton and curvaton fluctua- 
tions [77]. Since we assume the conventional quadratic 
curvaton potential, V(a) oc ra 2 <r 2 , we consider the latter 
possibility for which 



riNG 



1 



1 



1 + A- 



(26-7/ + ^), (53) 



See [69] for the full trispectrum including isocurvature contribu- 
tions. 
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where rj a cx H aa /H <C 77, since H 2 ^> p a during infla- 
tion. The projected la error on ting from future CMB 
missions is [78] 



An NG ~ 0.1 



50 

f local 
J NL 



(54) 



for a full-sky survey with Planck, and a factor of two 
smaller with CMBPol. The prospects of measuring 
a scale dependence are exciting because reconstruction 
may proceed without input from the amplitude /]v£ al , 
since ting by itself determines f 2 (cr). We consider such 
a detection in §111 E. 

Lastly, we discuss the possibility that neither isocur- 
vature modes nor non-Gaussianities are detected; in this 
case a precision measurement of the tensor spectral index 
can be leveraged to assist in the reconstruction effort. 



D. Tensor Index 

The energy density of the curvaton field during infla- 
tion is too small to generate a detectable tensor per- 
turbation, and hence the primordial gravitational wave 
spectrum is instead determined by the dynamics of the 
inflaton field, 



p h 



tit 



2 H 2 

IT 2 M|, ' 

-2e. 



(55) 
(56) 



However, since the curvaton does contribute to the over- 
all density perturbation, the tensor/scalar ratio is given 
by Eq. (22) and the single field consistency relation, 
r = — 8nT, is modified: 



-16tit 



2-f 2 (a)n T 



(57) 



A detection of r and tit can be used to break degen- 
eracies with single field inflation if their relationship is 
sufficiently different from the single field consistency re- 
lation r = — Stit [79, 80]. However, even if a particular 
observation confidently rules out single field inflation, it 
might not be clear what the alternative theory is since, in 
addition to curvatons, multifield [81, 82] and DBI infla- 
tion [17, 18], and even trans-Planckian effects [83, 84] lead 
to an altered consistency relation. An accurate identifi- 
cation of the underlying theory might require additional 
corroborating observations, like the amplitude and shape 
of non-Gaussianities or isocurvature modes. Our intent 
is not to enumerate all the different theories that predict 
the same relationship between r and ny, and so we will 
not view a measurement of tit as a degeneracy-breaking 
observation. Its status as such is mostly irrelevant for 
our purposes: as long as the measurements of r and tlt 
are consistent with the consistency relation in question 
to within experimental error, it can be used to constrain 



Reconstruction is possible because a determination 
of f(a) enables the inversion - four functions (V(0o), 
V((/>o), y'(0o), in terms of four observables 



CPs(fco), 



and tit). In fact, from Eq. (57) and 



Eqs. (23-25) we find that V /V cx y/e ex yJ—riT for 
the curvaton. Meanwhile, the single field reconstruction 
V /V cx y/e can be written either in terms of r or tlt 
by virtue of the single field consistency relation. The 
relative error between the curvaton and single field re- 
constructions of V /V is then 



A(V'/V) C 



Am 



A(V'/V) sing]e Ar 



(58) 



where the assumption is that r is better constrained than 
tit\ otherwise e would be most accurately determined by 
tit in both models, and the reconstructions would be 
equivalent. The other reconstruction diagnostic - the zo- 
ology classification - depends on the uncertainty in f 2 (a) : 



A/V) 



16 



+ 



r 1 + Ar/r nx 1 — Aut/tit 



(59) 



We see that when A(VyV) curv = A(V/ V) sing i e , the zo- 
ology reduces to that of single field inflation. An accu- 
rate detection of tlt will, however, be difficult with cur- 
rent technology: the Planck satellite in combination with 
current ground-based experiments like QUIET [85], BI- 
CEP [86], and PolarBear [87], might achieve a 1-a error 
of Atit ~ 0.1 [88]. While future space-based platforms 
such as CMBPol [89-91] might reduce this error by a 
factor of 2 in the absence of foregrounds [92] , a precision 
measurement of tit will most likely require the direct de- 
tection of primordial gravitational waves. 

The prospect of a direct detection of primordial grav- 
itational waves on scales ~ 0.1 — 1 Hz is the focus of 
several recent concept studies using space-based laser in- 
terferometers. Two proposals that have seen much ana- 
lytical attention are the Big Bang Observer (BBO) [93] 
and Japan's DECIGO project [94]. The forecasts estab- 
lished in these studies indicate that while BBO will be 
competitive with Planck (Atit ~ 0.1), upgrades to BBO 
might achieve Atit ~ 10 -2 , and DECIGO might reach 
Atit ~ 10 -3 [95, 96] in the most optimistic cases. These 
missions are still in the concept stage, and even if they 
become a reality, data will not become available for at 
least another decade. Even so, we include the possibility 
of a tensor index measurement for completeness. We now 
investigate the accuracy of the potential reconstruction 
that follows from a detection of each of the observables 
discussed in 5IIIB-IIID. 



E. Monte Carlo Analysis 

The mapping from observables to potential parame- 
ters, for example Eqs. (15 - 17) in the case of slow roll 
inflation, is a lowest order result. We wish to obtain a 
higher-order reconstruction; however, while the task of 
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inverting the observables (r, n s , dn s /dlnk, • • • ) to obtain 
the flow parameters (e, 77, A2, • • • ) in order to reconstruct 
the potential (via Eqs. (6-8)) is tractable at low order, 
it becomes prohibitively difficult at higher order. This 
difficulty motivated the development of Monte Carlo re- 
construction [97] in which one begins with the flow pa- 
rameters: e, 77, A2, • • • , Am, taken as initial conditions, 
and for each randomly drawn set obtains the inflationary 
evolution by solving the set of flow equations [98, 99], 



dH 

dN 
de 
dN 
d\j 
dN 
d\nk 
~dN 



eH, 

2e(77 -e), 

[(^-I)r7-^]A, + A m , *G[1,M], 

6-1, 



(60) 



where the system is truncated by taking Am+i = 0. 
Our time variable, dN = —Hdt, is the number of efolds 
before the end of inflation. In traditional applications 
[31, 33, 39, 98, 99], only solutions that yield sufficient 
inflation to solve the horizon and flatness problems are 
retained, typically with |A7V| G [45,70]. However, in 
this study, we impose the relatively diminished restric- 
tion |A7V| > 10, as this corresponds to the time period 
during which length scales directly probed by cosmolog- 
ical experiments exit the horizon. We work to 6 t/l -order 
in slow roll, and we have verified that our conclusions are 
robust up to 10 £/l -order. For each inflationary solution, 
observables and potential coefficients are separately de- 
termined. There is no need for a difficult inversion; the 
flow parameters are stochastically sampled and models 
corresponding to an observation of interest can then be 
selected out. The result is a consistent, high-order recon- 
struction that improves upon the analytic reconstruction 
given by Eqs. (15-17). This approach, of course, gener- 
alizes to include the curvaton. In order to obtain a ro- 
bust determination of possible degeneracies, it is impor- 
tant to sample a wide array of different inflation models 
and compare their potential reconstructions under the as- 
sumptions of both single field inflation and the curvaton 
model. 

We perform two flow analyses: one for single field in- 
flation in which the inflaton is soley responsible for gen- 
erating primordial perturbations, and the other for the 
curvaton scenario, in which the curvaton contributes to 
primordial perturbations. For the case of single field in- 
flation, we draw the initial flow parameters uniformly 
from the ranges 



[0,0.8], 
10^ +1 [ 



-0.5,0.5], 



(61) 



and we take lnfe^ = —8.047 so that the largest scale cor- 
responds to the quadrupole. We solve the flow equations 
and if |A7V| > 10, we calculate the potential parameters 



Eqs. (6-8), and the observables [100], 

r = 16e[l -2C(7/-e)], 
n s - 1 = 2r] - 4e - 2(1 + C)e 2 



(62) 



dn s 
dink 



-1(3-5^)677 + ±(3-C)A 2 , (63) 
4e(r/ - e) - 4e 2 - 12er/ + 2A 2 , (64) 



n T = -2c - (3 + C)e 2 + (1 + C)er/, 



(65) 



at k = O.OlMpc" 1 . Here C = 4(ln2 + 7 )-5 « 0.0815 and 
7 is the Euler-Mascheroni constant. We have included 
the running of the spectral index in order to obtain a 
higher-order reconstruction. 

For the curvaton scenario, we draw the flow parameters 
from the ranges Eq. (61) and uniformly sample f(cr) G 
[25, 10,000]. The lower bound on f(cr) is chosen so that 
the limit Eq. (31) is accurate, which can be verified by 
consulting the numerical solution for f(<j) in [36], while 
the upper bound is arbitrary. We solve the flow equations 
and if |A7V| > 10, we calculate the potential coefficients 
Eqs. (23-25) and the observables, 



16e[l - 2C(rj - e) 
2r/ - 2e 



-2e + 



l + /»6 



2(1 + C)e 2 



-i(3-5C)677 +i(3-C)A 2 , 



(66) 



(67) 



dn s 4e 2 - 12er] + 2A 2 4/ 2 (a)e(77 - e) 2 
— 4e(77 - e) ^7— h 



dink 



1 + P(a)e 



(l + / 2 (a)e) 2 ' 
(68) 

n T = -2e - (3 + C)e 2 + (1 + C)er], (69) 

at k = O.OlMpc" 1 . We now present results for the dif- 
ferent observational outcomes and collect our findings in 
Table 1. 



1. No Detection of Non-Gaussianity or Isocurvature Modes 

We first make projections for reconstruction utiliz- 
ing expectations from both the Planck Surveyor and 
CMBPol for the case in which neither isocurvature per- 
turbations nor non-Gaussianities are observed. We do 
not include the possibility of a precision measurement 
of tit at this time. This can be considered the worst- 
case scenario in which the degeneracy persists, and shall 
serve as a baseline against which to compare different 
observational outcomes. For Planck, we assume 68% 
CL detections of r (r > 0.01, Ar - 0.03) [101], n s 
(An s - 0.0038), and dn s /d\nk (Adn s /d\nk - 0.005) 
[102]. For CMBPol, we assume 68% CL detections of 
r (r > 10~ 4 , Ar - r/10), n s (An s - 0.0016), and 
dn s /dlnk (Adn s /dlnk - 0.0036) [90]. We will assume 
this base set of observations throughout the remainder 
of this analysis, unless otherwise indicated. The tensor 
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(d) 



(e) 



(f) 



FIG. 2. Monte Carlo results of single field (magenta) and curvaton (black) reconstructions in the absence of non-Gaussianities, 
isocurvature modes, or a precision measurement of tit- We present results for three fiducial values of r: r — 0.005 (a),(d), 
r = 0.05 (b),(e), r = 0.15 (c),(f). Top row presents results for Planck and bottom row for CMBPol. 



spectral index will not be adequately constrained with 
these CMB missions and the modified consistency rela- 
tion will not be useful for constraining curvaton models. 



We performed a 6^-order flow analysis on both single 
field inflation and the curvaton scenario, collecting 50,000 
models of each. We plot our results in Figure 2 for three 
fiducial values of r: r = 0.005 (a,d), r = 0.05 (b,e) and 
r = 0.15 (c,f), with Planck on top and CMBPol on bot- 
tom. Since (a) lies below the Planck detection threshold, 
in this case Planck can only impose the upper bound 
r < 0.01. The magenta points represent single field mod- 
els and the black points curvaton models. Figure 2 is 
a robust, higher-order representation of the degeneracies 
first mentioned in the context of Eqs. (23-25); the uncer- 
tainty in the potential parameters in the curvaton sce- 
nario is a result of our inability to constrain f 2 (cr). One 
can view the curvaton as 'contaminating' the single field 
result - by allowing for the presence of curvatons in the 
reconstruction, constraints are degraded relative to single 
field inflation by a factor of around four for V /V and a 
factor of five for V" /V. We assume fiducial values for the 
spectral index and running: n s =0.97 and dn s /dlnk = 0, 
although there is little sensitivity to the values of n s and 
dn s /d\nk chosen: changing n s acts to shift the distribu- 
tions slightly up and down, while changing dn s /d\nk has 



no noticeable effect. 9 

The inclusion of the curvaton does not affect the re- 
construction of the height of the potential, V(0o)- The 
higher-order reconstruction is consistent with the con- 
clusion drawn analytically at lowest-order. The lack of a 
tensor detection by Planck (r < 0.01) will therefore only 
impose an upper bound onV(^o), but does not otherwise 
strongly affect the errors on V /V and V" /V , as seen by 
comparison of Figure 2 (a) with (b). If Planck fails to 
constrain r, future concepts like CMBPol offer the tan- 
talizing possibility of detecting B-modes with r > 10 -4 . 
In such a case (r = 0.005), the single field reconstruction 
is notably improved (magenta models in Figure 2 (d)): 
V f /V, which is directly proportional to y/r, improves by 
a factor of 5 over the case in which Planck fails to detect 
r. However, surprisingly, such a detection has virtually 
no effect on the curvaton reconstruction (black models in 
Figure 2 (d).) This is generally true at the other fiducial 
r values as seen by comparison of Figures 2 (b) with (e) 
and (c) with (f), due to the poor constraints on f 2 (cr) 



9 This is because V ' (4>o) and V"((j)o) depend only weakly on run- 
ning (V n/ ((f)o) depends more strongly); it is therefore unlikely 
that even an accurate determination of dn s /d\nk : as might be 
possible with future 21-cm line observations [103], will be help- 
ful. 
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in the degenerate case; the percentage error on V /V for 
curvatons is 



v/v 



oc 



2l Ar\\(AP(a) 



r J 



1/2 



(70) 



While Ar/r < 1 for CMBPol, the error on f 2 (cr) is large: 
Af 2 (cr)/ j(<j) ~ 4000, and easily overwhelms the con- 
straints on r. We thus reach an important conclusion: in 
the degenerate case, the sole benefit of a detection of r 
will be the determination of the energy scale of inflation; 
the constraints on V /V and V" /V are largely insensi- 
tive to a detection or lack of detection of r, even by more 
advanced probes like CMBPol. 10 Nevertheless, in what 
follows we will continue to include r as a base observ- 
able because it will become important when we consider 
measurements of nr §111 E 4. 

As a result of the curvaton degeneracy, a given obser- 
vation is consistent with a greater variety of potentials. 
Recall that in single field inflation, for example, mod- 
els that predict r < — S(n s — l)/3 satisfy V" < and 
(logV)" < and are classified as small field potentials. 
In light of the effect of the degeneracy, it is now un- 
clear whether regions in the observable n s -r plane map 
uniquely to classes of functions of similar form. The ex- 
tent to which an unresolved curvaton preserves this map- 
ping and the resulting model classification provides an 
alternative view of the degeneracy problem, complemen- 
tary to the constraints found on V((j>) in Figure 2. We 
generate a new zoology in the presence of curvatons by 
assigning each of the models in Figure 2 to a zoology 
class (small field, large field, hybrid) as per the discus- 
sion at the end of §11. We first apply this methodology 
to single field models, and reproduce in Figure 3 (a) the 
zoology introduced in Figure 1. The slight overlap in re- 
gions is a result of the higher-order reconstruction of our 
analysis. In Figure 3 (b), we present the modified zo- 
ology that results when the curvaton is also present. In 
both cases we collect models within the observable ranges 
set by the WMAP7 95% CL: r < 0.55, n s G [0.91,1.2], 
and dn s /d\nk G [-0.09,0.02] [46]. The zoology is most 
clearly presented with each class plotted separately, so 
as to avoid overlap. We find that a substantial region of 
parameter space is degenerate in classification. In par- 
ticular, all observables compatible with large field models 
are also found to be consistent with hybrid models. The 
boundary between small and large field models lies at 
V /7 (0o) = 0- At lowest order, this occurs when e = —r] 
(c.f. Eqs. (6-8)), which gives a boundary in the observ- 



able plane: 
r = 12/~ 2 - 2r ± 2f~ 2 



1/2 



(71) 



where r = n s — 1. The boundary evidently depends on 
the value of /: as / — >> oo, this curve becomes more hor- 
izontal, approaching the line r = 0. While only a lowest 
order approximation, this result gives some intuition as 
to why the large field models extend all the way to the 
line r — in Figure 3 (b). 

The same approach applied to the boundary between 
hybrid and large field models, given by (logV^))" = 0, 
and thus e = r/, leads to the relation 



16r 



pr-2 



(72) 



This curve also tends to the line r = as / — » oo, and 
coincides with the single field case in the limit / — » 0. 
This agrees with our second finding, that all observations 
compatible with small field models are compatible with 
both large field and hybrid models. Hybrid models popu- 
late the full observable parameter space, extending down 
through the large field and small field regions. 11 Only 
those hybrid models existing in the single field 'hybrid' 
region can be correctly classified in the presence of the 
curvaton, i.e. they must satisfy r > 8(1 — n s ). For per- 
spective, we present the latest marginalized constraints 
on n s and r from WMAP+BAO+# [46] in Figure 4. 
Only observables outside the yellow region, a minority 
of the high CL parameter space, can be uniquely clas- 
sified according to the inflationary zoology. Note that, 
as in single field models, a detection of observables with 
r > 8(1 — n s )/3 would rule out small field inflation mod- 
els. 

It is particularly notable that large field curvaton mod- 
els can accommodate a small tensor amplitude on ob- 
servable scales. Recall that in single field inflation, the 
relation 



d(j) 
dN 



M P1 V2e, 



(73) 



together with r = 16e imposes a lower bound on the 
variation of the inflaton field during inflation, 



Ac/) > Mpi a - A TV, 



(74) 



where AN corresponds to observable scales. A tensor 
amplitude satisfying r < 0.01 indicates that A(j) < Mp\ 
while observable scales left the horizon. 12 In the presence 



In this statement we are assuming that background pollution of 
the tensor signal can be accurately subtracted. Detection of r is, 
of course, significant in that it provides unprecedented evidence 
for the existence of gravitational waves. 



11 In [104] an explicit hybrid potential is studied that yields 
r = 0, n s < 0, typical small field observables, using curva- 
tons/modulated reheating. 

12 Of course, it is still possible that A<t> > Mpi over the full course 
of inflation, even when r < 0.01 [105, 106]. 
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(a) (b) 



FIG. 3. (a) Zoology of single field models in the n s -r plane, obtained using the flow formalism (c.f. Figure 1). Blue, green, and 
red models represent hybrid, large field, and small field potentials, as defined in §11. (b) Zoology of models in which a curvaton 
is present but remains undetected. Because the classes overlap, we plot each model-type separately. Small field models are 
completely degenerate with large field models. Both types are completely degenerate with hybrid models. 



of the curvaton the bound is modified: 

As f 2 (cr)r — > 16, the field variation grows without bound. 
From Eq. (22), this correspond to the limit of small r: 
in the presence of curvatons, r < 0.01 no longer implies 
that Acf) < Mp\ even across observable scales. 13 



2. Detection of CDM Isocurvature Modes 

We now extend the results of the previous section by 
considering a detection of an isocurvature contribution 
to the density perturbation. We assume the same de- 
tections of r, n S: and dn s /dlnk : and analyze the three 
cases discussed in §111 B: an uncorrelated CDM mode 
(/? = 0), an anti-correlated CDM mode (/? = —1), and 
a CDM mode with an arbitrary correlation angle. We 
present results of the Monte Carlo in Figure 5 for the 
same fiducial values of r as considered previously. We 
plot the new constraints in orange and we include the 
constraints that are obtained in the absence of an isocur- 
vature detection (Figure 2) for comparison. It should be 
clear that the detection of an isocurvature mode rules out 
single field inflation and breaks the degeneracy; however, 
we will routinely make use of the constraints obtainable 



Information about the field range may be recovered if f(cr) can 
be accurately constrained. 
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FIG. 4. WMAP7+BAO+# 68% and 95% CL marginalized 
contours in the n s -r plane together with estimated Planck 
error bars. The (red) line segments show the predictions for 
V((f)) — m 2 (j) 2 and V(cj)) — A(/> 4 for the number N of e- folds 
before the end of inflation at which a mode crossed outside 
the horizon in the range N — [50, 60]. Only models predicting 
observables outside the yellow region (black solid line) can be 
uniquely classified according to the inflationary zoology when 
an unresolved curvaton is present. Small field models lie below 
the dashed black line. 



on single field inflation in the absence of any degeneracy- 
breaking observables as a benchmark against which to 
gauge the quality of the reconstruction. 

The most well-constrained case is that of the uncor- 
related CDM mode, Figure 5 (a) and (b). In fact, this 
case is as well-constrained as single field inflation; the 




FIG. 5. Monte Carlo results (orange) incorporating a detection of CDM isocurvature perturbations: uncorrelated (a) and (b), 
arbitrary correlation (c) and (d), anti-correlated (e). Models from Figure 2 §111 El, in which an isocurvature contribution is 
absent, are included for comparison. 
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uncertainty in the amplitude |<Scdm/C| 2 m (^2) a ^ 
Planck precision (~ 4% [47]) is subdominant relative to 
the errors on the other power spectrum observables. The 
function f(a) is well-constrained: Af 2 (cr) ~ 5. 

We present the results of the arbitrarily correlated 
CDM mode in Figure 5 (c) and (d). We have chosen 
f3 = —0.6 as a fiducial value with a Planck-precision error 
Af3 ~ 0.04 [47]. The constraints on V((j>) remain signifi- 
cantly improved relative to the case in which isocurvature 
modes are not detected, but deteriorate slightly relative 
to a detection of uncorrelated modes, with A/ 2 (a) ~ 250. 

The results of a detection of an anti-correlated mode 
are presented in Figure 5 (e). There is only one obser- 
vational window to portray, since r — >> in the limit 
that the isocurvature mode becomes anti-correlated (c.f. 
Eq. (45)). In this limit, f 2 (cr)e — >> oo, and the poten- 
tial reduces to the form Eqs. (28-30). Since the func- 
tion f(cr) drops out of the reconstruction, the errors on 
V /V and V" /V are due only to the uncertainties in the 
spectral parameters r, n s , and dn s /dlnk. The curvaton 
reconstruction is consequently comparable in precision 
to single field inflation in the absence of an isocurvature 
detection. Lastly, since the detection of an isocurvature 
component rules out single field inflation, there is no de- 
generacy problem, and no need to revisit the zoology 
classification in this case. 



3. Detection of N on- Gaussianity 

We next extend the results of §111 E 1 by considering 
a detection of non-Gaussianities. We again assume the 
same detections of r, n s , and dn s /dlnk as previously; we 
do not consider a detection of isocurvature modes in this 
section. Even with the most optimistic limits on the non- 
linearity parameters: Af l ° c L al - 1 (68% CL) [62, 63] and 
Ar^ al - 560 (95% CL) [71], we find that a detection 
of both the bispectrum and trispectrum does not improve 
constraints relative to those found in §1/7 E 1 for the case 
in which non-Gaussianities are not detected. The rea- 
son is simply that the projected uncertainty in rj^ 1 is 
too large to sufficiently constrain r a , with the result that 
even a precision detection of f 1 ^ 1 is incapable of suffi- 
ciently constraining / 2 (cr), with a resulting uncertainty 
Af 2 (cr) « 1500. This is an example of a case where the 
degeneracy is broken but the inversion problem persists. 

While measurements of f^L and tnl do not enable a 
successful inversion, a measurement of the scale depen- 
dence of non-Gaussianities by itself can be utilized in the 
reconstruction. Since Aung ~ l//iVL (Eq. (54)), the 
most optimistic outcome is when f 1 ^ 1 saturates its cur- 
rent observational bound, /jyff 1 ~ 70. We consider this 
case with the four fiducial values: ung = —0.2, —0.1, 0, 
and 0.1 at Planck precision (A%g = 0.07). We present 
results in Figures 6 (a) and (b) and note the improve- 
ment over the degenerate case Figure 2 (b) and (c) . For 
smaller /jy^ 1 , the uncertainty in ung grows and this ob- 



servable ceases to be of any help when f 1 ^ 1 < 35. A 
detection with CMBPol will further improve on these re- 
sults, with a minimum 

/tvl 11 > 1 7 - while the constraints 
on V(<fr) improve in these cases, we find that f 2 (cr) is not 
sufficiently constrained to enable a unique model classi- 
fication according to the zoology. 



4- Measurement of nr 

The last case that we consider is an extension of §111 E 1 
to include a precision measurement of the tensor spectral 
index, which will impose constraints on f 2 (cr) through 
the consistency relation Eq. (57). We consider separately 
Planck and CMBPol detections of r, and investigate the 
improvement in reconstruction that results from a di- 
rect detection of primordial gravitational waves, as might 
be possible with future space-based laser interferometers. 
The accuracy with which these probes will determine ut 
on direct detection scales oc 0.1 — 1 Hz is [95, 96] 



6 x 10- 18 
XQ GW h 2 ' 



(76) 



where VLqw is the gravity wave density, h is the present 
day Hubble parameter (in units of 100 km s _1 Mpc -1 ) 
and X characterizes the particular experiment: X = 0.25 
(BBO-standard), X = 2.5 (BBO-grand), and X = 100 
(DECIGO). The quantity flcwh 2 can be related to the 
primordial tensor spectrum [107], 



^Gwh 2 = A GW P h (k), 



(77) 



where the transfer function, Aqw, which governs the evo- 
lution of the perturbations up through horizon re-entry, 
has the numerical value Aqw = 2.74 x 10 -6 . We take 
the tensor spectrum to be of the form 



Ph(k) = P h (k ) 



(78) 



where olt = dnr/dlnk is the tensor index running. Al- 
though expected to be small and unlikely to be con- 
strained in future experiments, the running becomes im- 
portant when considering the large disparity between the 
scales of CMB probes (fco = 0-01 Mpc -1 ) and laser inter- 
ferometers (here taken to be k* = 6.5 x 10 14 Mpc -1 ). 
Our intent is to utilize this measurement of ut on direct 
detection scales to constrain the value of nr at the pivot 
where the consistency relations are valid. However, one 
cannot simply extrapolate the error Eq. (76) to CMB 
scales since this presumes that n^^o) = ^t(^*) 5 which 
does not hold in the presence of running. The uncertainty 
in ut at the pivot is therefore 



Am 



6 x 10- 18 
XA GW P h (h)_ 



-a T ln 



1/2 



(79) 
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r = 0.05 r = 0.15 




0.05 0.1 0.15 0.2 0.25 0.05 0.1 0.15 0.2 0.25 

V'M p /V VM p /V 

(a) (b) 



FIG. 6. Monte Carlo results for the case in which a scale dependent non-Gaussianity is detected in the CMB at Planck precision. 
For each fiducial value of r we consider four fiducial dependencies: ung — —0.2 (purple), —0.1 (yellow), (orange), and 0.1 
(red). From Figures 2 (b) and (c), we include the degenerate curvaton (grayed-out for clarity) and single field models (magenta) 
for reference. 



where olt — 4e?7 — 8e 2 . Since the running is un- 
known, it acts as an additional source of uncertainty 
when trying to match n^(^) to Ut^q). The run- 
ning is not, however, completely free to vary. It must 
be such that the tensor amplitude on direct detec- 
tion scales gives the correct amplitude of gravitational 
waves, £Igw, to within experimental error of the fidu- 
cial value. The la error on VLgw is approximately set 
by the signal-to- noise ratio of the experiment, AVLqw ~ 
(S/N)~\ with S/N = Xn GW /10- 18 [95]. We take 
as our fiducial values nx(&*) = —r/S and Vtcwh 2 — 
AGwrP®(ko)(k*/ko) nT ( k *\ consistent with a power law 
tensor spectrum. 

As discussed previously, as long as r and tit are consis- 
tent with Eq. (57) to within experimental error, the rela- 
tion can be used in the reconstruction. For ease of com- 
parison, we choose to base our reconstructions on fiducial 
values for r and tlt that agree with the single field re- 
lation; other choices will result in similar constraints on 
/ 2 (<j), since for a given value of r, Aut depends only 
weakly on the actual value of tit chosen. 

We present Monte Carlo results for Planck in Figure 
7: (a) and (b) correspond to fiducial values (r = 0.05, 
n T = -0.00625) and (r = 0.15, n T = -0.01875), respec- 
tively, and we compare the reconstruction projections of 
a direct detection of tensors with BBO-standard, BBO- 
grand, and DECIGO. We retain the curvaton (black) 
and single field (magenta) models from Figures 2 (b) 
and (c) and overlay the reconstructions from each di- 
rect detection experiment. For r = 0.05, direct detec- 
tion offers no improvement for single field inflation since 
Ar < Anx- However, there is improvement in V jV for 
curvatons: BBO-standard (orange), BBO-grand (red), 
and DECIGO (purple). Overall, constraints are im- 



proved for r = 0.15, where we find that single field un- 
certainties are reduced with DECIGO (cyan). The im- 
provement for curvatons is similar to that seen in Figure 
7 (a), and the curvaton and single field reconstructions 
are almost equivalent with DECIGO. For both curvatons 
and single field, we note the rather insignificant improve- 
ment of DECIGO over BBO-grand, despite the factor of 
forty increase in sensitivity offered by DECIGO. We con- 
sider the projections for CMBPol in Figure 7 (c)-(e). At 
r = 0.15, the errors on the curvaton reconstruction are 
reduced to those of single field inflation (purple points 
behind cyan). Constraints are weakest for r = 0.005 - 
BBO-standard offers no improvement over the degener- 
ate case. 

It is important to emphasize that the improvement 
in reconstruction that results from the direct detection 
of gravity waves relies on optimistic outcomes of mis- 
sions still in the planning stages, and should therefore be 
viewed as a best-case scenario. In the future, the ten- 
sor index might instead be constrained by a more accu- 
rate measurement of B-modes on CMB scales. For com- 
pleteness, we consider projections of an ideal detection of 
B-modes on CMB scales in Figures 7 (f)-(h). The con- 
straints are competitive with direct detection, achieving 
an accuracy close to that of CMBPol-j-BBO-grand. The 
uncertainties in tit for the ideal experiment are taken to 
be An T = 0.01 at r = 0.005, An T = 0.009 at r = 0.05, 
and An T = 0.007 at r = 0.015 [108]. 

Given the marked improvement in constraints that a 
direct detection affords, the expectation is that some of 
the uniqueness of model classification, characteristic of 
single field inflation, might be retained even in the pres- 
ence of curvatons. The uncertainty in the function f 2 (cr) 
determines an approximate zoology through Eqs. (71) 
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FIG. 7. Monte Carlo results for single field and curvaton models with a direct detection of tensors at three different precisions: 
BBO-standard, BBO-grand, and DECIGO, for Planck (top row) and CMBPol (middle row). We present results at three fiducial 
values: (r = 0.005, n T = -0.00065), (r = 0.05, n T = -0.00625), and (r = 0.15, n T = -0.01875). We include the curvaton 
and single field models without a measurement of nr and overlay the reconstructions from direct detection at each different 
precision. For curvatons: BBO-standard (orange), BBO-grand (red), DECIGO (purple); for singe field we only present results 
from DECIGO (cyan), since the other reconstructions do not improve on the case without direct detection (magenta). In the 
last row we present projections in which ut is determined by an ideal detection of B-modes on CMB scales rather than a direct 
detection experiment: curvatons (green) and single field (blue). 



and (72), and an accurate measurement of ut can be 
used to constrain it through the consistency relation Eq. 
(57). In Figure 8, we present the zoology for the cases 
considered in Figure 7. Regions in which the degener- 
acy persists - where different classes of single field and 
curvaton models make the same observable predictions - 
are shaded in gray. With Planck, the variation in f 2 (cr) 
is large for r = 0.05 for each tit detection, and the zo- 
ology is degenerate (Figure 3 (b)) with hybrid models 
the only unique class. While the degeneracy is reduced 



with Planck+DECIGO at r = 0.15 (Figure 8 (a)), there 
is little room for a unique classification: the la error 
on r spans the entire y-axis, so small field models can- 
not be uniquely classified, while large field models might 
be identified for a precise measurement of n s « 0.955. 
With CMBPol, however, the zoology for r = 0.05 is 
improved with BBO-grand and DECIGO (Figure 8 (b) 
and (c)); with DECIGO, there exist regions in which 
a unique classification of hybrid, small, and large field 
models is possible. CMBPol at r = 0.15 is shown in 
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Planck w/ DECIGO, r = 0.15 



' 0.93 0.94 0.95 0.96 0.97 0.98 0.9' 



(a) 



CMBPol w/ BBO-grand, r = 0.05 



CMBPol w/ DECIGO, r = 0.05 



0.93 0.94 0.95 0.96 0.97 0.98 0.99 1 



0.93 0.94 0.95 0.96 0.97 0.98 0.99 



(b) 



(c) 



CMBPol w/ BBO-standard, r = 0.15 



CMBPol w/ BBO-grand, r = 0.15 



CMBPol w/ DECIGO, r = 0.15 






(d) 



(e) 



(f) 



FIG. 8. Zoology in the presence of the curvaton degeneracy for different quality measurements of r and 71t, at r = 0.05 and 
r = 0.15. Degenerate regions are shown in gray. Experiments not shown are fully degenerate. 



Figures 8 (d)-(f). Unsurprisingly, the best combination 
is CMBPol+DECIGO for which the zoology almost re- 
duces to that of single field inflation. Lastly, we note that 
for r = 0.005, the zoology is degenerate for all direct de- 
tection sensitivities as well as an ideal B-mode detection. 



IV. DISCUSSION 

We have investigated the degeneracy problem that 
arises when the primordial power spectrum is generated 
by something other than the inflaton. We analyzed the 
curvaton as an example of this class of models. The cur- 
vaton was selected for our case study because it has been 
analyzed extensively in the literature and is phenomeno- 
logically similar to other scenarios in this class, such as 
modulated reheating. While the curvaton is representa- 



tive of these models, we have analyzed but one corner of 
the potentially degenerate model space, and our results 
should therefore be considered a best-case scenario for re- 
construction in the absence of any degeneracy-breaking 
observable. The importance of further input from the- 
oretical model building to help diminish the degeneracy 
and inversion problems should not be underestimated: 
the more information we have about the underlying fun- 
damental physics behind inflation, the better off we are 
in terms of the reconstruction program. For example, if 
the inflaton has DBI Lagrangian and is non-minimally 
coupled to gravity, interesting observational signatures 
can arise [109]. Under certain circumstances the degener- 
acy between minimally and non-minimally coupled DBI 
inflationary models can be resolved [110]. We consider 
another major source of possible degeneracies that can 
arise from non-canonical inflation in [19]. 
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TABLE I. Reconstruction results of each case considered in this analysis. In the first set of columns the errors on the potential 
coefficients in the presence of curvatons are given relative to those expected from single field inflation; in cases where the 
observation rules out single field inflation, the values are given for reference only. In the second set of columns these errors are 
given relative to the worst-case degeneracy, i.e. no detections of non-Gaussianity, isocurvature modes, or nr. When constraints 
depend on the fiducial value of r chosen, we provide results separated by a slash, with the constraints for r = 0.005 followed 
by those for r — 0.05 and r = 0.15 (for Planck only these latter two values are relevant); if only one number is provided then 
there is no difference. In the last column, 'All' indicates that all three zoology classes can be uniquely reconstructed; otherwise, 
only those classes that can be uniquely reconstructed are listed. NG: non-Gaussianities, iso: isocurvature, H: hybrid, LF: large 
field, SF: small field. 



Here, we first determined the size of the degeneracy by 
performing Monte Carlo reconstructions on both single 
field and curvaton models in the absence of observables 
that might serve to distinguish the curvaton, such as lo- 
cal non-Gaussianities and/or isocurvature modes. We 
chose .P$(fco)> r i n si an d dn s /d\nk as base observables 
common to both models and found that reconstruction 
is degraded in the presence of an unresolved curvaton 
because it precludes a unique inversion of the observ- 
ables to obtain the potential. We then introduced ad- 
ditional observables into our study: isocurvature modes, 
local non-Gaussianities, and the tensor spectral index. 
While the former two observations break the degeneracy 
by ruling out single field inflation and lending support 
to the curvaton, our main focus was to study how well 
a set of degeneracy-breaking observables can be inverted 
to reconstruct the scalar potential. 

A comprehensive view of our results is presented in 
Table 1. We list the size of the errors on V /V and V" /V 
that result from a detection of each observable considered 
in the analysis. The errors are given relative to the errors 
expected from single field inflation (in the event that the 
observation rules out single field inflation, constraints are 
compared to single field as a benchmark), and in the 
second set of columns we display the uncertainties from 
each observation relative to the totally degenerate case. 
The uncertainty in f 2 (cr) is most directly relevant to the 
zoology portraits shown throughout the paper. Our main 
conclusions are as follows: 

• Null-detection. In the absence of a detection of 



non-Gaussianities or isocurvature modes - the de- 
generate case - constraints on V /V and V" /V are 
each degraded by a factor of 4 and 6.5, respectively, 
relative to single field inflation. In terms of the zo- 
ology, we find that only models with observables 
satisfying r > —S(n s — 1) can be uniquely classi- 
fied. These are hybrid models (see Figure 4). Even 
the improved measurements of r and n s possible 
with CMBPol offer scant improvement. 

• Non-Gaussianity. The amplitude of non- 
Gaussianities, /jyjf 1 , serves as a useful way to dis- 
criminate between the two theories, but the es- 
timator is otherwise unhelpful for reconstruction. 
Even with single field inflation ruled out in this 
case, a detection of more than the amplitude of 
the bispectrum is needed for inversion. We in- 
cluded the possibility of a measurement of the 
trispectrum but found that the projected uncer- 
tainties of future missions are too large to im- 
prove constraints on V(<j)) beyond the degenerate 
case. However, a measurement of the scale depen- 
dence of non-Gaussianities offers modest improve- 
ment with Planck when /J°£ al > 35 (/#£ al > 17 
with CMBPol.) 

• Isocurvature. A detection of an isocurvature contri- 
bution to the primordial density perturbation also 
serves to break the degeneracy and the amplitude 
and correlation angle of the mode facilitate a suc- 
cessful reconstruction. For a CDM isocurvature 
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component, an uncorrected mode offers the best 
constraints. The zoology becomes virtually identi- 
cal to that of single field inflation. 

• Gravitational waves. A precise measurement of the 
tensor index, tit, enables the use of the consistency 
relations of the degenerate models. Even when the 
values of r and tlt satisfy both the single field and 
curvaton relations and the degeneracy persists, a 
precision measurement of tit can substantially re- 
duce the degeneracy. The detection must be of suf- 
ficient quality: An T < G>(10~ 2 ), with r > 0.01 of- 
fering the best reconstructions; this precision might 
be achieved with a future direct detection of pri- 
mordial gravitational waves. Since the curvaton re- 
construction makes use of the consistency relation, 
a quality measurement of r is also advantageous. 
With a Planck detection of r, the degeneracy be- 
gins to clear for fiducial r = 0.15 with DECIGO, 
but the degeneracy remains for smaller r. How- 
ever, with a CMBPol detection of r, the degen- 
eracy resolves partially for r > 0.01 with BBO- 
grand and almost completely with DECIGO. In 
the case of small tensor/scalar ratio, r < 0.01, 
we find that with CMBPol+BBO-grand and DE- 
CIGO, constraints on V(<fr) improve over the case of 
no tit but the zoology remains degenerate. Lastly, 
an ideal detection of B-modes on CMB scales of- 
fers a reconstruction capability similar to that of 
CMBPol+BBO-grand. 

In summary, a detection of an isocurvature component, 
the scale dependence of local non-Gaussianities, or a pre- 
cision measurement of the tensor spectral index enable 
successful reconstructions. The consistency relation is 
most powerful when the constraints on r and tlt are opti- 
mized - with CMBPol and DECIGO, respectively. While 
the reconstruction in this case is practically equivalent 
to single field inflation, we require a favorable observa- 
tional outcome (r > 0.01), and optimistic foregrounds. 



We end by considering the possibility that r < 0.01 and 
goes unmeasured by Planck. In the degenerate case, a 
CMBPol detection of r < 0.01 will constrain the energy 
scale of inflation, but it will not, by itself, improve con- 
straints on V /V or V" /V. In this case, the prospects 
of experimentally probing the detailed physics of the in- 
flationary mechanism are rather dire. While detection 
of r will provide enthralling evidence for the existence of 
gravitational waves, we have discovered that in order for 
CMBPol to adequately reconstruct the shape of V ((/>), it 
must be paired with a precision direct detection of pri- 
mordial gravitational waves. 

The resolution of the degeneracy and inversion prob- 
lems will require a herculean effort on the part of preci- 
sion observational cosmology: we will need better probes 
of B-mode polarization, advanced direct probes of pri- 
mordial gravitational waves, and an improved under- 
standing of foregrounds. It is important that we under- 
stand what we can expect to constrain with these future 
missions. In this analysis we have taken a small step 
towards attaining this understanding. 
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